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ABSTRACT 

Near-Infrared  spectroscopy  combined  with  high  spatial  resolution  imaging  have  been  used  in  this 
work  to  probe  the  central  500  pc  of  M82.  Imaging  observations  in  the  2.36  jim  CO  band  head  are  added 
to  our  previously  published  near-infrared  hydrogen  recombination  line  imaging,  near-infrared  broadband 
imaging,  and  3.29  jim  dust  feature  imaging  observations,  in  order  to  study  the  nature  of  the  starburst 
stellar  population.  A  starburst  model  is  constructed  and  compared  with  the  observations  of  the  stellar 
clusters  in  the  starburst  complex.  Our  analysis  implies  that  the  typical  age  for  the  starburst  clusters  is 
107  yr.  In  addition,  our  high  spatial  resolution  observations  indicate  that  there  is  an  age  dispersion 
within  the  starburst  complex  that  is  correlated  with  projected  distance  from  the  center  of  the  galaxy.  The 
inferred  age  dispersion  is  6  x  106  yr.  If  the  starburst  in  M82  is  propagating  outward  from  the  center,  this 
age  dispersion  corresponds  to  a  velocity  of  propagation,  originating  in  the  center,  of  ~  50  km  s_1.  Our 
quantitative  analysis  also  reveals  that  a  Salpeter  initial  mass  function,  extending  from  0.1  to  100  M0, 
can  fit  the  observed  properties  of  M82  without  using  up  more  than  30%  of  the  total  dynamical  mass  in 
the  starburst. 

Subject  headings:  galaxies  evolution  —  galaxies:  individual  (M82)  —  galaxies  nuclei  — 
galaxies:  starburst  —  galaxies:  star  clusters  —  infrared:  galaxies 


1.  INTRODUCTION 

Numerous  galaxies  have  been  observed  to  be  undergoing 
an  intense  period  of  star  formation  activity  that  completely 
dominates  the  output  radiation  from  the  preexisting  galaxy. 
Although  there  is  general  agreement  on  the  source  for  the 
observed  energetic  activity  in  many  of  these  galaxies,  funda¬ 
mental  questions  on  the  causes  of  starburst  episodes,  their 
progression  and  evolution,  and  the  detailed  properties  of 
the  stars  formed  in  them,  are  not  definitely  answered. 
Efforts  to  gain  a  firm  understanding  of  some  of  these 
properties  will  be  aided  greatly  by  detailed  high  spatial 
resolution  observations  of  nearby  starburst  galaxies. 

M82,  close  by  and  hence  well  studied,  has  served  as  the 
archetypical  starburst  galaxy,  providing  an  ideal  testing 
ground  in  which  to  study  starburst  theories.  In  this  paper, 
starburst  models  are  constructed  and  compared  with  a 
series  of  detailed  high  spatial  resolution  near-infrared 
imaging  observations  described  in  depth  by  Satyapal  et  al. 
(1995),  hereafter  Paper  I,  and  Greenhouse  et  al.  (1997).  In 
Paper  I,  high  spatial  resolution  near-infrared  Fabry-Perot 
imaging  observations  of  the  central  kiloparsec  of  M82,  in 
conjunction  with  near-infrared  broadband  imaging  obser¬ 
vations  and  imaging  in  the  3.29  jim  dust  feature,  were  used 
to  examine  the  extinction  toward  the  starburst  region,  the 
state  of  the  ionized  gas,  and  the  nature  of  the  stellar  popu¬ 
lation. 
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In  this  article,  we  conduct  a  quantitative  analysis  of  the 
major  points  discussed  in  Paper  I.  In  addition,  new  imaging 
observations  of  the  2.3  jim  CO  absorption  feature  are  added 
to  our  previously  published  results,  providing  further 
insight  into  the  starburst  episode  in  M82.  Unlike  most  past 
studies  of  M82,  we  model  the  properties  of  the  individual 
stellar  clusters  in  the  central  500  pc  region. 

In  the  next  section,  we  describe  our  new  CO  band  head 
observations  of  M82,  and  the  results  are  presented  in  §  3. 
The  stellar  clusters  within  the  central  500  pc  are  identified, 
and  their  major  intrinsic  properties  are  summarized  in  §  4. 
Our  starburst  model  is  described  in  §  5  and  is  applied  to  the 
stellar  clusters  in  M82  in  §  6. 

2.  CO  2.3  jim  ABSORPTION  FEATURE  OBSERVATIONS 
2.1.  Background 

Observations  of  the  first-overtone  bands  of  CO  near  2.3 
jim  have  been  used  widely  to  probe  the  stellar  population  in 
galaxies.  This  feature,  observed  only  from  the  atmospheres 
of  cool  stars,  shows  a  systematic  increase  in  depth  with 
increasing  stellar  luminosity  and  decreasing  effective  tem¬ 
perature  (Frogel  et  al.  1978,  hereafter  F78;  Kleinmann  & 
Hall  1986).  It  is  also  found  to  increase  with  increasing 
metallicity  (Frogel,  Persson,  &  Cohen  1983).  The  strength  of 
the  feature  is  defined  by  the  “CO  index,”  which  is  essen¬ 
tially  the  difference  in  magnitudes  between  the  center  of  the 
CO  obsorption  and  the  extrapolated  continuum. 

Observers  have  used  a  variety  of  spectral  resolutions  and 
wavelengths  to  measure  the  CO  index.  In  the  early  work  on 
CO  indices  by  F78,  the  stellar  CO  absorption  feature  was 
characterized  using  the  “photometric  CO  index,”  COF, 
defined  as  the  difference  in  magnitudes  between  the  narrow- 
band  filters  centered  at  2.2  jim  (Al  =  0.11  jim)  and  2.36  jim 
(A2  =  0.08  jim)  relative  to  the  star  a  Lyrae.  Unlike  these 
measurements,  our  present  observations  (A2  =  0.04  jim) 
permit  the  measurement  of  the  slope  of  the  continuum, 
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Fig.  1  c  Fig.  Id 

Fig.  1. — 2.36  gm  CO  absorption  distribution  compared  with  the  Bry  emission  in  M82.  (a)  CO  index  image  with  contours  overlaid.  Contours  displayed 
from  0.1  to  0.3  arcsec-2  in  intervals  of  0.05  arcsec-2.  Tick  marks  correspond  to  5".  North  is  up,  and  east  is  to  the  left.  (6)  IF(Bry)  image  with  contours  overlaid. 
Contours  displayed  from  1.5  to  2.0  arcsec-2  in  intervals  of  0.1  arcsec-2.  (c)  Bry  image  with  contours  overlaid.  Contours  displayed  from  1.5  x  10-17  to 
3  x  10- 17  Wm- 2  arcsec-2  in  intervals  of  5  x  10-18  Wm-2  arcsec-2  (d)  JC-band  image. 


allowing  any  necessary  correction  for  deviations  from  a  flat 
continuum,  e.g.,  reddening.  In  addition,  the  higher 
resolution  observations  allow  us  to  choose  wavelengths 
that  avoid  contamination  of  the  measured  flux  from  strong 
emission  lines  that  may  be  present  in  the  spectrum  of 
objects  such  as  M82  (e.g.,  Bry ;  l  =  2.166  /tm). 

Observations  of  the  CO  index  in  the  dust-enshrouded 
nuclei  of  starburst  galaxies  can  serve  as  a  powerful  tool  in 
the  study  of  the  late-type  stellar  population.  It  has  the 
advantage  of  being,  in  principle,  independent  of  extinction, 
mass,  and  distance  and  can  thus  provide  potentially  an 
observationally  rigid  constraint  to  models  of  the  starburst 
population. 

2.2.  CO  Band  Head  Images  of  M82 :  Observations  and  Data 
Reduction 

The  CO  band  head  observations  of  M82  were  carried  out 
in  1995  June  at  the  Wyoming  Infrared  Observatory 
(WIRO)  with  the  university  of  Rochester  Third  Generation 
Array  Camera  using  the  1.85%  bandwidth  cold  circular 
variable  filter  (CVF).  Images  were  taken  with  the  CYF  set 
to  2.140,  2.250,  and  2.360  jim.  Ten  20  s  on-source  integra¬ 
tions  and  10  background  integrations  at  each  of  these  wave¬ 
length  settings  were  obtained.  On-chip  nodding  was  used  to 
double  our  efficiency.  The  plate  scale  was  0''24  pixel-1,  and 
the  spatial  resolution  was  1'.'2.  The  wavelengths  were  chosen 
after  examining  a  high-resolution  spectrum  of  an  M5  III 


star  from  the  Kleinmann  &  Hall  (1986)  stellar  atlas  and 
optimizing  on  the  depth  of  the  absorption  feature  for  our 
spectral  resolution.  The  continuum  settings  were  selected  to 
measure  the  slope  accurately  and  to  ensure  that  these  wave¬ 
lengths  were  not  near  emission  lines  that  would  be  present 
in  M82.  A  standard  star  with  no  CO  absorption  in  its  spec¬ 
trum  was  chosen  as  a  flux  calibrator  (v  UMa;  F2  IV,  mK  = 
2.99).  The  two  continuum  images  were  used  to  determine 
the  extrapolated  continuum  at  2.36  /im,  and  an  image  of  the 
CO  index  was  generated.  Using  an  electronic  version  of  the 
Kleinmann  &  Hall  (1986)  stellar  atlas,  the  CO  index  was 
calculated  synthetically  using  our  spectral  resolution  and 
choice  of  wavelengths  and  using  those  of  F78,  and  the 
results  were  compared.  Our  present  CO  index  was  found  to 
be  roughly  consistent  with  J.  A.  Frogel’s;  the  relationship 
between  the  two  is  given  by 

COF  =  1.127  COthiswork-  0.004.  (1) 

3.  DISTRIBUTION  OF  THE  2.36  gm  CO  ABSORPTION 
FEATURE  IN  M82 

Normal  galaxies  have  typical  CO  indices  of  0.15,  with 
very  little  scatter  around  this  value  (F78).  Deeper  absorp¬ 
tion  depths  in  galaxies  can  be  produced  by  a  large  popu¬ 
lation  of  late-type  supergiants  and  massive  giants,  or  by  a 
metal-enriched,  otherwise  normal  population  of  stars.  Shal¬ 
lower  depths  can  be  produced  by  an  enhanced  population 
of  early-type  stars  or  by  the  contamination  of  the  K-band 
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continuum  from  nonstellar  processes  (e.g.,  free-free  and  free- 
bound  emission  from  ionized  gas,  and  the  emission  from 
hot  dust).  The  image  of  the  CO  index  of  M82,  together  with 
our  Bry  image  (Paper  I),  and  our  image  of  the  Bry  equiva¬ 
lent  width,  fY(Bry),  are  shown  in  Figure  1.  A  K-band  image 
is  also  displayed  for  comparison.  The  CO  index  ranges  from 
~  0.26  in  individual  pixels  to  ~  0  across  the  central  500  pc. 
As  can  be  seen  in  Figure  1,  the  CO  index  is  distinctly  anti¬ 
correlated  with  the  Bry  equivalent  width  in  M82;  the  CO 
index  is  high  at  the  nucleus  where  IF(Bry)  is  low  and  is  at  its 
minimum  ~  12"  west  and  3"  south  of  the  nucleus  where 
IF(Bry)  is  at  its  maximum.  Since  the  Bry  emission  is  tracing 
both  the  bluest  stellar  population  and  the  regions  where  the 
contribution  to  the  K-band  continuum  from  nebular  and 
dust  emission  is  highest,  this  anticorrelation  is  expected. 

4.  INTRINSIC  PROPERTIES  OF  THE  STELLAR  CLUSTERS 

In  this  section,  we  summarize  the  major  properties  of  the 
stellar  clusters  found  in  the  central  500  pc  of  M82.  In  order 
to  determine  the  intrinsic  properties  of  the  stellar  clusters  in 
the  M82  starburst,  an  extinction  correction  must  be  applied 
to  the  observed  emission.  In  Paper  I,  we  used  an  extinc¬ 
tion  map  that  was  generated  using  the  observed  Pa/f/Bry 
flux  ratio,  assuming  case  B  recombination,  and  the  t;  oc 
i-1-85  (Landini  et  al.  1984)  extinction  law  to  deredden  the 
broadband  images,  generating  extinction-corrected  J,  H, 
and  K  images.  We  also  calculated  the  free-free  and  free- 
bound  components  of  the  broadband  emission  in  an  effort 
to  separate  out  the  stellar  component. 

4.1.  Identification  of  Pointlike  Sources 

The  extinction-corrected  K-band  image  was  used  to 
search  for  pointlike  sources  using  the  task  DAOFIND  in 


TABLE  l 


Positions  of  the  12  Pointlike  Broadband  Near-Infrared 
Sources  Identified  in  the  Central  500  pc  of  M82 


Source  Number 

Right  Ascension  Offset” 

Declination  Offset” 

1 . 

-6.7 

2.6 

2 . 

-4.6 

1.4 

3 . 

-3.5 

-2.6 

4 . 

0.0 

0.0 

5 . 

3.8 

-2.6 

6 . 

5.0 

-0.1 

7 . 

7.1 

-0.6 

8 . 

0.4 

-3.9 

9 . 

0.9 

-1.7 

10 . 

10.8 

-3.5 

11  . 

10.9 

-1.8 

12 . 

13.2 

-3.4 

”  Offsets  from  K  nucleus  in  arcseconds.  Eight  of  the  sources  corre¬ 
spond  to  sources  identified  by  Pipher  et  al.  1987 


the  DAOPHOT  software  package  (Stetson  1987).  This 
algorithm  searches  for  local  maxima  in  the  input  image 
whose  amplitudes  are  greater  than  a  user-specified  detec¬ 
tion  threshold.  Twelve  pointlike  sources  were  identified  and 
are  displayed  in  Figure  2,  and  their  positions  are  given  in 
Table  1.  Based  on  their  luminosity  and  spatial  extent  (Table 
2)  at  the  distance  of  M82,  these  pointlike  sources  are  likely 
to  be  young  compact  star  clusters.  In  order  to  determine  the 
properties  of  these  clusters,  2"  diameter  circular  apertures 
were  centered  on  each  of  the  12  sources,  and  standard  aper¬ 
ture  photometry  was  carried  out  at  these  positions.  Initial 
attempts  to  extract  the  magnitudes  of  these  pointlike 
sources  involved  experimenting  with  PSF-fitting  algorithms 
in  the  DAOPHOT  software  package.  However,  these  algo¬ 
rithms  are  optimized  for  crowded  stellar  field  photometry 


Fig.  2. — Pointlike  sources  in  the  central  500  pc  of  M82.  Foreground  extinction-corrected  X-band  image  with  the  location  of  the  12  compact  objects. 
North  is  up,  and  east  is  to  the  left;  (0, 0)  corresponds  to  the  position  of  the  nucleus. 
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TABLE  2 

Summary  of  Properties  of  Pointlike  Sources  in  M82 


Source 

mJncor 

(•f-Jfluncor 

(H~K)aa  cor 

(H-K)  cor 

^ff  +  fb 
(%) 

^ff  +  fb 
(%) 

^ff  +  fb 
(%) 

*^bb 

(%) 

Kbb 

(%) 

^■dust 

(%) 

1 . 

10.8“ 

10.1“ 

1.03b 

0.77b 

0.1 7b 

0.20b 

1° 

2C 

3° 

id 

ld 

13e 

2 . 

10.6 

10.0 

1.02 

0.72 

0.28 

0.21 

2 

2 

3 

l 

1 

12 

3 . 

11.3 

11.0 

0.89 

0.55 

0.47 

0.28 

2 

2 

3 

l 

1 

13 

4 . 

9.9 

9.3 

1.40 

0.91 

0.53 

0.37 

1 

1 

1 

l 

0 

20 

5 . 

11.1 

10.1 

1.35 

0.94 

0.30 

0.27 

4 

4 

6 

2 

1 

15 

6 . 

11.1 

10.2 

1.31 

0.87 

0.28 

0.24 

3 

3 

5 

2 

1 

13 

7 . 

11.0 

10.2 

1.04 

0.71 

0.14 

0.13 

3 

3 

6 

2 

1 

7 

8 . 

11.4 

10.7 

1.23 

0.82 

0.41 

0.28 

4 

4 

6 

3 

1 

13 

9 . 

10.7 

10.0 

1.10 

0.80 

0.31 

0.31 

2 

2 

4 

1 

1 

19 

10 . 

11.5 

10.7 

1.23 

1.01 

0.27 

0.39 

7 

7 

11 

4 

2 

24 

11 . 

11.0 

10.3 

1.14 

0.80 

0.37 

0.33 

3 

3 

5 

2 

1 

18 

12 . 

12.0 

11.2 

1.26 

1.03 

0.35 

0.44 

8 

7 

11 

5 

2 

26 

a  We  estimate  a  10%  error  in  our  photometry.  All  quantities  correspond  to  a  2"  aperture  centered  on  specified  source. 

b  The  estimated  error  in  our  colors  is  15%.  The  corrected  colors  correspond  to  extinction-corrected  values.  The  free-free  and  free-bound 
components  were  subtracted  out. 

c  The  percent  contribution  from  free-free  and  free-bound  processes  at  the  indicated  bandpass. 
d  The  percent  contribution  from  the  Pa/1  and  Bry  line  flux  to  the  broadband  continuum. 

e  The  fraction  of  the  JC-band  continuum  from  thermal  emission  from  600  K  dust  required  to  generate  “  normal  ”  colors  (see  Fig.  3). 


with  minimal  overlying  nebulosity.  When  using  these  algo¬ 
rithms  on  the  extinction-corrected  K-band  image  of  M82, 
the  combination  of  luminous  and  variable  nebulosity, 
crowded  field  photometry,  and  the  lack  of  isolated,  pure 
point  sources  needed  to  characterize  the  point-spread  func¬ 
tion  (PSF)  resulted  in  unsuccessful  point-source  subtrac¬ 
tions.  We  therefore  used  standard  aperture  photometry  on 
our  images,  with  apertures  that  were  small  enough  to 
prevent  overlap  among  the  pointlike  sources  (see  Fig.  2). 
The  flux  in  these  apertures  was  taken  to  represent  the  total 
flux  from  the  stars  in  the  clusters,  the  emission  from  the 
ionized  gas,  and  dust  emission.  In  reality,  there  can  be  some 
small  contribution  to  the  continuum  emission  from  the 
underlying  old  stellar  population;  a  component  of  the  total 
ionized  gas  could  be  ionized  from  the  background  inter¬ 
stellar  radiation  field  and  not  entirely  by  the  stars  from  that 
cluster,  and  there  could  be  a  contamination  to  the  total 
emission  from  one  source  from  its  neighboring  sources. 


4.3.  Summary  of  Intrinsic  Properties  of  Pointlike  Sources 

Using  a  foreground  screen  extinction  model,  the  observed 
magnitudes  and  colors  of  each  of  the  point  sources  were 
corrected  for  extinction,  and  the  components  to  the  near- 
infrared  emission  described  in  Paper  I  were  calculated.  The 
results  are  summarized  in  Table  2.  In  addition  to  the  free- 
free  and  free-bound  components  to  the  broadband  contin¬ 
uum,  the  contributions  to  the  broadband  emission  from  the 
Pa/i  and  Bry  line  flux  were  also  calculated.  The  spectrum  of 
M82  includes  many  other  emission  lines,  but  they  do  not 
comprise  a  significant  component  of  the  broadband  emis¬ 
sion;  the  maximum  contribution  to  the  J-band  continuum 
from  the  Pa /?  line  flux  is  5%. 

Table  3  summarizes  various  other  properties  of  the  12 
pointlike  sources.  The  extinction-corrected  Bry  line  flux  was 
used  to  calculate  the  number  of  04  V  stars  (Panagia  1973) 
required  to  produce  this  flux;  the  emission  from  these  stars 


TABLE  3 

Summary  of  Properties  of  Pointlike  Sources  in  M82  Continued 


Source 

^Bry 

(1(T17  Wm'2) 

JVuv 

(s-1) 

^04V 

log  [W(Bry)  (A)] 

^"'^observed 

CO 

intrinsic 

1 . 

7.9 a 

6.1 

71b 

1.295c 

0.172d 

0.177e 

2 . 

10.5 

8.0 

94 

1.357 

0.206 

0.214 

3 . 

3.7 

2.8 

33 

1.311 

0.142 

0.147 

4 . 

8.1 

6.2 

73 

0.977 

0.215 

0.219 

5 . 

17.0 

13.1 

153 

1.637 

0.113 

0.121 

6 . 

12.9 

9.9 

117 

1.532 

0.208 

0.221 

7 . 

15.1 

11.6 

136 

1.597 

0.243 

0.261 

8 . 

10.6 

8.1 

96 

1.636 

0.154 

0.165 

9 . 

10.7 

8.2 

96 

1.363 

0.191 

0.199 

10 . 

18.4 

14.1 

166 

1.882 

0.129 

0.146 

11 . 

11.6 

8.9 

105 

1.536 

0.044 

0.047 

12 . 

11.2 

8.6 

101 

1.863 

0.157 

0.177 

a  We  estimate  a  10%  error  in  our  photometry.  All  quantities  correspond  to  a  2"  aperture  centered 
on  specified  source. 

b  The  number  of  04  V  stars  was  calculated  using  the  ionizing  photon  rates  from  Panagia  1973. 

c  ±0.01. 

d  The  estimated  uncertainty  is  0.015. 

e  The  intrinsic  CO  index  was  calculated  by  diluting  the  feature  depth  by  the  free-free  and  free- 
bound  processes  listed  in  Table  2. 
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in  the  K  band  was  found  to  contribute  less  than  1%  to  the 
total  K-band  continuum  at  each  of  the  sources.  The 
observed  CO  indices  of  each  of  the  sources,  along  with  the 
values  corrected  for  the  free-free  and  free-bound  continuum, 
are  also  listed  in  Table  3. 

The  observed  and  corrected  colors  of  the  12  sources  are 
shown  in  Figure  3.  Vectors  indicating  the  effects  of  fore¬ 
ground  extinction,  extinction  from  a  uniform  mixture  of 
dust  and  stars,  and  the  emission  from  600  K  dust  are  also 
shown.  The  colors  have  been  corrected  for  foreground 
extinction,  and  the  free-free  and  free-bound  components 
were  removed  (see  Paper  I).  As  can  be  seen  from  Figure  3, 
the  sources  have  corrected  colors  that  are  significantly  dif¬ 
ferent  from  those  of  normal  stars  or  clusters  of  normal  stars. 
All  have  H—K  colors  that  are  “too  red”  relative  to  their 
corresponding  J  — if  color.  Many  starburst  galaxies  have 
colors  that  occupy  this  region  of  the  color-color  diagram 
(see  Doyon,  Joseph,  &  Wright  1994;  Turner  1995;  Satyapal 
1995;  Smith  et  al.  1996).  If  we  assume  that  the  position  of 
these  sources  on  the  two-color  diagram  can  be  ascribed  to 
the  emission  from  a  normal  population  of  stars  and  the 
emission  from  dust,  we  can  calculate  the  fraction  of  the 
broadband  emission  attributable  to  dust  emission.  This 
fraction  is  listed  in  Table  2.  The  dust  emission  was  modeled 
by  graybody  dust  grains  with  emissivity  e  oc  2  2  (see 
Aaronson  1977)  and  a  temperature  of  600  K.  However, 
under  this  assumption,  the  dust  dilutes  the  CO  index,  and 
the  resulting  relatively  blue  colors  of  the  stellar  clusters  are 
generally  inconsistent  with  their  CO  indices  (see  Table  2)  for 
any  normal  population  of  stars.  The  CO  index  correspond¬ 


ing  to  the  J—H  colors  of  normal  supergiant  stars  is  dis¬ 
played  on  the  right  axis  in  Figure  3.  Thus,  we  conclude  that 
in  order  for  the  corrected  colors  to  be  consistent  with  the 
observed  CO  indices,  a  combination  of  purely  foreground 
extinction  and  dust  emission  cannot  be  the  only  effects  con¬ 
sidered. 

4.4.  Effects  of  Different  Models  of  the  Dust  Geometry 

In  the  previous  section,  we  have  shown  that  a  purely 
foreground  extinction  screen  does  not  result  in  extinction- 
corrected  colors  that  are  consistent  with  normal  stellar 
populations  and  that  dust  emission  cannot  be  invoked.  In 
this  section,  we  investigate  the  effects  of  varying  the 
assumed  dust  geometry  on  the  extinction-corrected  colors 
(assuming  no  dust  emission).  For  the  purpose  of  illustrating 
the  importance  of  these  effects  on  near-infrared  colors  in 
general,  we  will  use  the  integrated  colors  from  a  large  aper¬ 
ture  centered  on  the  M82  nucleus  as  representative. 

In  this  section,  we  will  consider  three  different  dust 
geometries,  all  considered  previously  in  the  study  of  M82 
(e.g.,  McLeod  et  al.  1993;  Puxley  1991;  Telesco  et  al.  1991). 
The  first  is  the  simple  foreground  screen  geometry  discussed 
above.  Here  the  ionized  gas,  toward  which  the  extinction  is 
calculated,  and  the  stars  are  being  obscured  by  a  single 
overlying  screen.  The  second  is  a  homogeneous  mixture  of 
stars,  gas,  and  dust.  Finally,  we  will  consider  a  two-screen 
geometry,  where  the  ionized  gas  and  stars  are  being 
obscured  by  a  dust  screen,  followed  by  another  emission 
region  of  stars  and  gas  emitting  a  fraction  of  the  first  emit¬ 
ting  region’s  intensity,  followed  by  another  dust  screen. 
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Fig.  3. — Two-color  diagram  of  the  M82  pointlike  sources.  The  filled  triangles  correspond  to  the  observed  colors,  while  the  filled  circles  correspond  to  the 
corrected  colors.  The  colors  have  been  corrected  using  the  foreground  screen-derived  recombination  line  extinction  map,  and  the  free-free  and  free-bound 
components  to  the  broadband  emission  have  been  removed.  The  three  dashed  lines  correspond  to  the  effects  of  foreground  extinction,  a  homogeneous 
mixture  extinction  geometry,  and  the  thermal  emission  from  600  K  dust.  The  tick  marks  on  the  reddening  vectors  correspond  to  xK  of  0.2.  The  tick  marks  on 
the  dust  emission  vector  correspond  to  5%  fractions  of  the  total  K  flux  density.  The  colors  of  normal  stars  are  also  shown  in  the  figure  (Koornneef  1983).  The 
CO  index  of  normal  supergiants  is  indicated  on  the  right  axis  (F78). 
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(II)  Homogeneous  Mixture  Geometry 
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These  geometries  are  illustrated  in  Figure  4. 7  The  resulting 
colors  are  shown  on  a  color-color  diagram  in  Figure  5.  We 
have  used  the  integrated  emission  from  a  24"  aperture  cen¬ 
tered  on  the  nucleus.  The  extinction  was  then  calculated  for 
each  of  the  discussed  dust  geometries  using  the  total  Bry 
and  Pa/1  line  fluxes  in  this  aperture  ( FBiy  =  4.24  x  10” 15  W 
m-2  and  FP.l/;  =  9.35  x  10“ 15  W  m-2),  adopting  the  same 
assumptions  of  case  B  recombination  discussed  earlier 
(Paper  I).  An  extinction  power  law,  rA  oc  2-1'85,  was  again 
employed.  The  resulting  values  for  tk  and  the  extinction- 
corrected  K-band  magnitudes  are  also  shown  in  Figure  5. 
In  addition,  we  have  calculated  the  CO  index,  also  dis¬ 
played  in  Figure  5,  from  the  integrated  flux  densities  at  2.36, 
2.25,  and  2.14  gn i  in  a  24"  aperture.  As  can  be  seen  from 


7  Adapted  from  McLeod  et  al.  (1993). 


Figure  5,  the  assumed  dust  geometry  has  a  significant  effect 
on  the  corrected  colors.  A  foreground  screen  geometry  does 
not  generate  corrected  colors  that  are  consistent  with  the 
colors  of  any  normal  stars.  However,  a  first-order  correc¬ 
tion  to  the  foreground  screen  geometry,  namely,  the  two- 
screen  geometry  with  /  =  0.3  and  g  =  0.2,  can  generate 
corrected  colors  characteristic  of  normal  stars  that  are  con¬ 
sistent  with  the  observed  CO  index  of  0.18.  Although  these 
two  models  generate  significantly  different  colors,  the 
extinction-corrected  K-band  magnitudes  are  not  appre¬ 
ciably  different.  The  difference  between  the  corrected  mag¬ 
nitudes  in  these  two  cases  is  only  0.24,  or  a  factor  of  1.2  in 
the  corrected  K-band  luminosities.  Thus,  a  small  variation 
of  the  foreground  screen  geometry  can  have  a  substantial 
effect  on  the  extinction-corrected  colors  of  the  M82  star- 
burst  without  significantly  affecting  the  corrected  magni¬ 
tudes.  Note  that  the  extinction  model  must  result  in 
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Fig.  5. — Two-color  diagram  of  the  representative  colors  of  the  M82 
starburst.  The  effects  of  different  dust  geometries  on  the  corrected  colors 
are  indicated.  Values  for  zK  and  mK  (before  and  after  extinction  corrections 
have  been  applied)  are  indicated.  The  average  value  for  the  CO  index  is 
also  shown. 


“normal”  colors  without  changing  substantially  the 
derived  corrected  magnitudes  from  their  corresponding 
foreground  screen  values.  As  mentioned  in  Paper  I,  the 
extinction-corrected  Bry  flux  cannot  overpredict  the 
thermal  bremsstrahlung  continuum  at  millimeter  wave¬ 
lengths.  Using  our  foreground  screen  extinction-corrected 
Bry  image,  the  total  extinction-corrected  flux  of  1.0  x  10“ 14 
W  m-2  results  (assuming  case  B  recombination,  Te  =  5000 
K)  in  a  derived  3  mm  free-free  continuum  flux  density  that  is 
already  ~75%  of  that  observed  (Carlstrom  &  Kronberg 
1991).  Since  the  emission  at  3  mm  has  some  nonthermal 
component  as  well  as  a  component  from  dust,  the 
extinction-corrected  Bry  flux  cannot  be  much  larger  than 
that  obtained  using  our  foreground  screen  extinction 
model.  Thus,  although  we  cannot  constrain  precisely  the 
extinction  geometry,  we  can  say  that  the  extinction  model 
must  result  in  “normal”  colors  and,  at  the  same  time, 
cannot  produce  corrected  magnitudes  substantially  differ¬ 
ent  from  those  obtained  using  a  foreground  screen  extinc¬ 
tion  model.  With  only  two  recombination  line  maps, 
however,  we  cannot  constrain  the  extinction  geometry  more 
precisely.  In  order  to  obtain  more  meaningful  extinction- 
corrected  colors  for  the  stellar  clusters  in  the  M82  starburst, 
at  least  one  additional  recombination  line  image  would  be 
required.  For  example,  imaging  in  the  Bra  (4.05  /mi)  line 
would  allow  us  to  conduct  a  more  rigorous  study  of  the 
extinction  in  the  M82  starburst,  resulting  in  colors  that 
would  be  more  representative  of  the  true  colors  of  the 
central  stellar  clusters. 

The  analysis  above  emphasizes  the  need  to  exercise 
caution  in  interpreting  extinction-corrected  colors.  For  the 
case  of  M82,  we  have  confirmed  a  predominantly  foreground 
screen  model  for  the  extinction,  but  we  have  shown  that  a 
slight  deviation  from  this  extinction  model  can  affect  the 
corrected  colors  significantly.  The  corrected  magnitudes, 
however,  are  not  affected  substantially.  Hence,  all 
extinction-corrected  absolute  quantities  are  well  deter¬ 
mined  in  our  study  of  M82.  The  extinction-corrected  colors, 
however,  are  not  determined  precisely  enough  for  compari¬ 


son  with  stellar  population  models,  for  which  very  small 
color  changes  are  significant. 

5.  STARBURST  MODEL 

In  order  to  conduct  a  quantitative  investigation  of  the 
properties  of  the  starburst  clusters  in  M82,  we  have  con¬ 
structed  our  own  set  of  population  synthesis  models.  Unlike 
most  previous  studies  of  starburst  galaxies,  our  models  can 
be  applied  not  only  to  the  integrated  properties  of  a  star- 
burst,  but,  in  the  case  of  M82,  to  the  observations  of  the 
individual  stellar  clusters. 

Our  model  is  analytic  and  assumes  a  power-law  initial 
mass  function  (IMF)  between  a  variable  lower  mass  cutoff 
and  a  fixed  upper  mass  cutoff  of  mv  =  100  M0.  The  other 
free  parameter  in  the  model  is  the  total  mass  assumed  in  the 
burst.  As  an  illustrative  case,  we  have  assumed  a  single  IMF 
slope  set  at  the  Salpeter  value  for  the  solar  neighborhood  of 
T  =  1.35  (Salpeter  1955).  More  recent  studies  of  the  solar 
neighborhood  IMF  have  shown  that  the  low-mass  IMF  is 
nearly  flat  (e.g.,  Bessell  &  Stringfellow  1993)  while  the  high- 
mass  IMF  for  field  stars  is  steeper  than  the  Salpeter  value 
(e.g.,  Massey  et  al.  1995b),  although  the  cluster  IMF  for 
high-mass  stars  in  the  solar  neighborhood  is  similar  to  the 
Salpeter  IMF  (Massey,  Johnson,  &  De  Gioia-Eastwood 
1995a).  We  have  taken  this  simple  case  to  illustrate  the 
relationship  between  the  observed  quantities  and  the 
parameters  derived  for  the  IMF.  The  stellar  tracks  for  solar 
metallicity  stars  (Z  =  0.02)  by  Schaller  et  al.  (1992)  were 
used.  We  have  also  assumed  an  instantaneous  burst  of  star 
formation.  In  reality,  the  star  formation  rate  in  a  starburst 
galaxy  nucleus  can  have  a  complicated  time  dependence. 
However,  unlike  previous  studies  of  M82,  we  are  modeling 
the  stellar  clusters  in  the  M82  starburst  that  are  character¬ 
ized  more  accurately  as  “instantaneous  starbursts”  than 
are  the  integrated  properties  from  the  entire  starburst 
region.  We  have  synthesized  the  J,  H,  and  A  magnitudes 
and  colors,  the  supernova  rate,  and  the  CO  index  for  a 
starburst  population  with  ages  ranging  from  105  to  108  yr. 
The  CO  index  was  calculated  using  the  analytical  expres¬ 
sions  from  Doyon  et  al.  (1994).  The  Bry  equivalent  width, 
JF(Bry),  from  Leitherer  &  Heckman  (1995)  was  compared 
with  our  data  also.  We  do  not  include  the  effects  of  a 
nebular  continuum  or  a  preexisting  population  of  stars.  Our 
data  permit  us  to  subtract  the  nebular  continuum,  allowing 
direct  comparison  with  the  model  quantities.  The  preexist¬ 
ing  stellar  population  does  not  contribute  significantly  to 
the  total  luminosity  from  the  central  30"  of  M82.  The  total 
foreground  extinction-corrected  A -band  magnitude  in  a  30" 
aperture  centered  on  the  nucleus  is  —22.0.  Assuming  that 
the  mass-to-K-band  luminosity  ratio  of  the  old  stellar 
population  is  comparable  to  the  value  of  23  Mo/L0 
observed  in  normal  galactic  bulges  (Devereux,  Becklin,  & 
Scoville  1987;  Thronson  &  Greenhouse  1988;  Gaffney, 
Lester,  &  Telesco  1993),  the  old  stars  would  have  an  absol¬ 
ute  A-band  magnitude  of  —18.1,  or  3%  of  the  total 
extinction-corrected  A-band  luminosity,  if  they  constituted 
50%  (3.5  x  10s  M0)  of  the  total  dynamical  mass  within  the 
central  450  pc.  The  results  of  our  model  were  found  to  be 
well  within  the  range  obtained  by  other  models  (e.g., 
Leitherer  &  Heckman  1995;  Chariot  &  Bruzual  1991). 

The  calculated  quantities  described  above  can  be  used  to 
constrain  the  age  and  IMF  of  a  starburst  when  compared 
with  observations.  We  now  compare  the  model  predictions 
with  the  starburst  episode  in  M82. 


No.  1,  1997 


PROPAGATING  STAR  FORMATION? 


155 


6.  APPLICATION  TO  M82 

6.1.  Age  of  the  Burst:  Propagating  Star  Formation! 

Several  observations  suggest  that  the  starburst  is  propa¬ 
gating  outward  from  the  center  of  M82.  The  A -band  emis¬ 
sion,  presumably  dominated  by  the  output  from  evolved 
stars,  is  relatively  smooth  and  centrally  concentrated  (see 
Fig.  2),  while  the  recombination  line  (Paper  I),  10  pm 
(Telesco  et  al.  1991),  and  3.29  pm  dust  feature  (Paper  I) 
emission  show  a  more  bilobal  distribution.  This  suggests 
that  the  nucleus  contains  an  older  stellar  population  and 
that  the  starburst  may  be  propagating  outward.  In  Paper  I, 
the  radial  variation  in  the  3.29  pm  dust  feature-to-  Bry  flux 
ratio  was  shown  to  provide  additional  qualitative  evidence 
for  the  presence  of  a  radially  dependent  age  variation  in  the 
central  500  pc.  Since  the  3.29  pm  dust  feature  is  thought  to 
be  the  result  of  nonequilibrium  heating  of  small  grains  by 
ionizing  and  lower  energy  photons,  while  the  Bry  emission 
is  related  only  to  ionizing  photons,  the  variation  of  the 
extinction-corrected  ratio  from  values  greater  than  35  near 
the  nucleus  to  6  at  the  position  15"  west  of  the  nucleus  may 
be  a  crude  indicator  of  the  hardness  of  the  radiation  field 
within  the  starburst  region,  again  suggesting  a  radially 
dependent  evolutionary  effect.  In  this  section,  a  more  quan¬ 
titative  investigation  of  the  age  of  the  starburst  complex 
using  our  starburst  models  will  be  carried  out. 

The  age  of  the  stellar  clusters  within  the  M82  starburst 
can  be  inferred  most  accurately  by  the  observed  CO  index 
and  IF(Bry),  as  has  been  done  previously  for  other  starburst 
galaxies  (e.g.,  Doyon  et  al.  1994).  These  observations  are 
particularly  well  suited  to  constraining  the  age  of  a  stellar 
population  since  the  model  quantities  show  a  sharp  change 
in  a  short  time  interval.  The  evolution  of  the  CO  index  and 
IF(Bry)  for  a  starburst  model  with  mL  =  1  M0  and  mv  = 
100  M0,  together  with  the  observed  values  at  the  positions 
of  the  12  pointlike  sources  discussed  in  §  4,  is  shown  in 
Figure  6.  These  model  quantities  are  relatively  insensitive  to 
the  lower  mass  cutoff  assumed  and  thus  can  be  used  to 
derive  the  age  of  the  starburst  independently.  The  observed 
values  have  not  been  corrected  for  dust  emission  or  the 
small  continuum  contribution  from  an  old  stellar  popu¬ 
lation.  These  components  would  not  have  a  substantial 
effect  on  the  CO  index  or  much  less  the  Bry  equivalent 
width.  The  CO  indices  imply  ages  between  4  x  106  and 


1  x  107  yr  for  the  12  different  sources,  roughly  consistent 
with  the  ages  inferred  from  the  observed  Bry  equivalent 
widths.  Although  the  precise  age  depends  on  the  detailed 
time  dependence  of  the  star  formation  rate,  the  observed 
CO  indices  and  Bry  equivalent  widths  imply  a  relatively 
narrow  range  in  age  (the  effects  of  these  model  quantities  on 
the  assumed  star  formation  rate  can  be  found  in  Doyon  et 
al.  1994).  Thus,  the  implied  age  for  the  M82  starburst  is 
approximately  107  yr.  The  most  recent  estimate  of  the  age 
of  the  starburst  in  M82  by  Rieke  et  al.  (1993)  is  (1.3- 
3)  x  107  yr,  depending  on  the  details  of  the  assumed  star 
formation  rate.  Their  model  of  the  integrated  emission  from 
the  central  ~  500  pc  required  incorporation  of  two  bursts. 
We  have  obtained  estimates  of  the  age  of  the  stellar  clusters 
directly.  At  these  times,  the  most  massive  stars  have  com¬ 
pleted  nuclear  burning,  and  the  effective  temperature  of  the 
ionizing  stars  is  ~  35,000  K,  consistent  with  estimates 
derived  from  various  forbidden  line  ratios  (e.g.,  Achtermann 
&  Lacy  1995).  Early  starburst  models  imposed  an  upper 
mass  cutoff  to  account  for  the  relatively  soft  radiation  field; 
since  the  effective  temperature  of  the  ionizing  stars 
decreases  naturally  with  time  as  the  most  massive  stars 
evolve  off  the  main  sequence,  this  assumption  is  not  neces¬ 
sary. 

Apart  from  determining  the  age  of  the  M82  starburst  as  a 
whole,  our  high  spatial  resolution  data,  together  with  our 
starburst  models,  allow  a  quantitative  investigation  of  the 
possible  radial  age  gradient.  In  this  scenario,  star  formation, 
which  initially  started  at  the  nucleus,  induces  star  formation 
at  other  locations  as  the  newly  formed  population  of  stars 
compress  the  surrounding  regions.  With  Figure  6,  we  can 
estimate  the  age  dispersion  implied  by  the  observed  W(Bry) 
and  CO  index  within  the  central  ~  500  pc.  The  estimated 
age  dispersion  obtained  through  the  spread  in  CO  indices  of 
the  12  different  sources  is  ~6  x  106  yr.  The  maximum  pro¬ 
jected  separation  is  ~  230  pc,  implying  a  velocity  of  propa¬ 
gation,  originating  in  the  center,  of  40  km  s-1.  Using  the 
spread  in  the  values  of  W{ Bry),  the  inferred  velocity  is  60  km 
s_1.  This  velocity  is  not  an  unreasonable  velocity  for  expan¬ 
sion.  Achtermann  &  Lacy  (1995)  find  evidence  for  radial 
expansion  of  the  ionized  gas  at  ~  55  km  s-1  from  the 
nucleus  of  M82.  The  correlation  between  the  CO  index  and 
W(Bry)  for  the  12  pointlike  sources  is  shown  in  Figure  7. 

The  dependence  of  the  CO  index  and  the  Bry  equivalent 


Fig.  6. — The  CO  index  as  a  function  of  time  for  a  starburst  model  with  mL  =  1  M0  and  mv  =  100  M0,  together  with  the  intersections  of  the  observed 
values  of  the  12  pointlike  sources. 
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Fig.  7. — CO  index  vs.  W( Bry)  for  the  12  pointlike  sources 

width  on  the  projected  radius  for  the  12  different  pointlike 
sources  is  shown  in  Figure  8.  There  is  a  general  anti¬ 
correlation  of  the  CO  index  with  the  Bry  equivalent  width 
as  a  function  of  projected  radius;  radii  close  to  the  nucleus 
are  generally  associated  with  larger  CO  indices  and  lower 
Bry  equivalent  widths,  while  the  opposite  situation  prevails 
at  larger  radii.  This  anticorrelation  is  most  pronounced 
along  the  major  axis  of  the  galaxy,  as  can  be  seen  in  Figure 
9.  It  should  be  noted  that  the  decrease  in  the  CO  index  with 
radius  is  not  the  result  of  the  “  filling  in  ”  of  the  absorption 
feature  with  increasing  nonstellar  continuum,  as  can  be  seen 
from  Figure  10.  In  addition  to  correcting  the  CO  index  for 
dilution  from  the  nebular  component  discussed  in  §  4.3  and 
Paper  I,  the  dilution  resulting  from  “hot”  dust  emission 
was  also  estimated.  The  near-infrared  emission  from  a  col¬ 
lection  of  small  particles  transiently  heated  by  the  absorp¬ 
tion  of  ultraviolet  and  visible  photons  resembles  that  of  a 
graybody  with  color  temperature  1000  K  (Sellgren,  Werner, 
&  Dinnerstein  1983).  Sellgren  et  al.  have  shown  that  the 
3.29  /im  line-to-continuum  ratio  is  roughly  constant  at 
several  locations  in  a  number  of  reflection  nebulae.  Using 
this  constant  ratio  of  ~  6,  the  extinction-corrected  3.29  /an 
feature  image,  and  the  assumption  that  the  hot  dust  contin¬ 


-100  -50  0  50  100 

Radial  Distance  (pc) 


Fig.  9. — CO  index  and  W(Bry)  along  the  major  axis  of  the  galaxy 


uum  can  be  approximated  by  a  1000  K  graybody,  the 
extinction-corrected  flux  density  from  hot  dust  emission 
near  2  /an  was  calculated  and  accounted  for  in  the  CO 
index  image.  Source  7  appears  to  have  an  anomalously  high 
CO  index  compared  with  its  Bry  equivalent  width.  We  note 
that  our  CO  band  head  and  Bry  data  were  obtained  using 
the  same  telescope  and  camera,  and  the  data  were  reduced 
consistently.8  Our  measured  CO  index  at  the  position  of  the 
nucleus  using  the  appropriate  apertures  was  found  to  be  in 
good  agreement  with  that  found  by  McLeod  et  al.  (1993) 
and  Lester  et  al.  (1990).  The  scatter  in  the  CO  index  versus 
the  fT(Bry)  graph  is  therefore  not  a  reflection  of  the  obser¬ 
vational  uncertainty  of  our  data.  The  high  CO  index  of 
source  7  may  be  due  to  a  higher  metallicity  in  this  source, 
although  the  required  metallicity  would  be  extreme.  Using 
the  calibration  derived  from  the  observations  of  globular 
clusters  (Frogel  et  al.  1983),  if  the  observed  CO  index  of  0.24 
in  source  7  is  the  result  of  a  metal-enriched  population  of 

8  Error  bars  for  these  quantities  are  displayed  in  Fig.  8  and  listed  in 
Table  3. 


Projected  Radius  (pc)  Projected  Radius  (pc) 


Fig.  8. — Left  panel :  CO  index  for  the  12  different  sources,  identified  by  number,  vs.  projected  radius  from  the  nucleus.  The  estimated  uncertainty  in  the 
CO  index  is  shown  in  the  lower  right-hand  corner.  Right  panel :  W(Bry)  vs.  projected  radius.  The  estimated  uncertainty  is  shown  in  the  lower  right-hand 
corner. 
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Fig.  10. — Radial  variation  of  the  observed  (filled  circles )  and  intrinsic 
(i open  circles)  CO  index  along  the  major  axis  of  M82.  The  intrinsic  CO 
index  was  calculated  by  removing  the  effects  of  the  free-free,  free-bound, 
and  hot  dust  emission  near  2  /an.  An  estimate  of  the  hot  dust  emission, 
calculated  by  using  the  3.29  /an  feature  flux,  and  the  assumption  of  a 
line-to-hot  dust  continuum  ratio  of  6  (Sellgren  et  al.  1983),  was  also  sub¬ 
tracted  out. 


otherwise  normal  giants,  the  metallicity  would  have  to  be  a 
factor  of  ~  12  larger  than  solar  (CO  =  0.16  +  0.074[Fe/H], 
where  [Fe/H]  is  the  logarithm  of  the  metal  abundance  rela¬ 
tive  to  solar;  Frogel  et  al.  1983).  The  metallicity  dependence 
of  the  CO  index  was  determined  from  data  on  globular 
cluster  giants  with  below  solar  abundances ;  the  reliability  of 
this  calibration  for  the  present  data  is  therefore  unknown. 
Although  the  correlation  is  not  remarkable,  these  data 
suggest,  most  importantly,  that  there  is  an  age  dispersion 
within  the  starburst  region  of  M82,  and  that  this  age  disper¬ 
sion  appears  to  be  correlated  generally  with  projected 
radius.  Most  previous  starburst  models  have  treated  the 
integrated  quantities  within  the  inner  ~450  pc  as  a  single 
burst.  Rieke  et  al.  (1993)  approximated  the  continuum  of 
starburst  activity  by  two  bursts  in  their  models  in  order  to 
account  for  roughly  two  characteristic  populations — a 
young  blue  one  farther  from  the  nucleus  and  an  older, 
redder  nuclear  population.  As  we  have  seen,  there  is  a  con¬ 
tinuum  of  starburst  activity  possibly  moving  outward, 
probably  as  a  result  of  sequential  star  formation. 

6.2.  Is  the  IMF  Deficient  in  Low-Mass  Stars  ? 

Our  extinction  corrections  (Paper  I)  yield  an  extinction- 
corrected  K-band  luminosity  in  a  30"  aperture  centered  on 
the  nucleus  of  M82  that  is  a  factor  of  at  least  1.6  smaller 
than  previous  values.  What  effect  does  this  have  on  the 
inferred  IMF? 


In  this  section,  we  will  treat  the  integrated  properties 
from  a  30"  (450  pc)  aperture  centered  on  the  nucleus  as  a 
single  starburst  population.  We  have  shown  that  there  is 
actually  an  age  dispersion  among  the  stellar  populations 
within  this  aperture.  It  is  therefore  a  crude  approximation 
to  treat  the  entire  starburst  region  as  a  single  instantaneous 
burst  characterized  by  a  single  age;  this  analysis  is  intended 
only  to  illustrate  the  relationship  between  the  observed 
quantities  and  the  parameters  derived  for  the  IMF,  with  an 
emphasis  on  testing  the  validity  of  a  low-mass  deficient  star- 
burst  population  in  M82.  These  characteristics  are  critically 
dependent  on  the  extinction  correction  applied  to  the 
K-band  luminosity  and  the  total  mass  assumed  in  the  burst. 
Neither  of  these  quantities  is  known  unequivocally.  In  order 
to  make  firm  conclusions  about  the  IMF  in  the  M82  star- 
burst,  we  stress  the  importance  of  interpreting  the  predic¬ 
tions  of  these  models  in  the  light  of  the  current 
observational  uncertainties. 

The  extinction-corrected  K-band  magnitude  of  M82  has 
varied  widely  in  the  literature.  These  variations  are  due  to 
differences  in  the  K-band  photometry  and  to  differences  in 
the  adopted  extinction  geometry.  Table  4  is  a  summary  of 
the  observed  and  extinction-corrected  (where  available)  K 
magnitudes  in  the  literature  from  a  30"  (450  pc)  aperture 
centered  on  the  nucleus.  We  have  demonstrated  in  Paper  I 
that  our  recombination  line-derived  extinction  is  the  most 
reliable  estimate  for  the  extinction  to  date  toward  the 
central  starburst  region  of  M82.  We  have  also  shown  (see 
§  4.4)  that  a  slight  modification  to  the  foreground  screen 
assumption,  a  modification  that  is  necessary  to  explain  the 
corrected  colors  of  the  starburst  region  in  terms  of  the  emis¬ 
sion  from  a  normal  population  of  stars,  does  not  change 
substantially  the  extinction-corrected  K-band  magnitude. 
Thus,  Mk  in  a  30"  aperture  centered  on  the  nucleus  of  M82 
is  very  close  to  our  value  of  —22.0  (Paper  I)  listed  in  Table 
4. 

The  total  mass  involved  in  the  burst  is  also  not  deter¬ 
mined  uniquely.  The  dynamical  mass  within  the  starburst 
region  consists  of  the  total  mass  of  the  gas,  the  total  mass  of 
stars  involved  in  the  burst,  and  the  mass  of  stars  from  an  old 
population  of  stars  predating  the  burst.  There  is  some 
observational  uncertainty  in  the  estimates  for  the  dynamical 
mass  within  the  starburst  region.  There  is  also  a  large 
uncertainty  associated  with  the  fraction  of  this  mass  that  is 
molecular,  and  finally  the  mass  of  the  preexisting  stellar 
population  is  not  known.  Hence,  there  is  some  flexibility  in 
the  mass  that  can  be  assigned  to  the  starburst  population. 

The  dynamical  mass  has  been  estimated  using  visible 
(e.g.,  Gotz  et  al.  1990)  and  mid-infrared  emission  lines  from 
the  ionized  gas  (e.g.,  Beck  et  al.  1978;  Achtermann  &  Lacy 
1995),  stellar  absorption  lines  (e.g.,  Gaffney  et  al.  1993),  and 
molecular  lines  (Lo  et  al.  1987,  Shen  &  Lo  1995).  These 
estimates  can  vary  by  as  much  as  a  factor  of  1.5.  Table  5 
summarizes  a  selection  of  these  observations.  In  order  to 


TABLE  4 

Selection  of  K-Band  Measurements 


Authors 

Observed 

Mk 

Extinction-corrected  MK 

Extinction  Model 

Rieke  et  al.  1980  . 

-23.3 

Homogeneous  mixture,  derived  from  H  recombination  lines 

Telesco  et  al.  1991 . 

-22.1 

-22.6 

Foreground,  derived  from  H  —  K  color 

McLeod  et  al.  1993 . 

-22.5 

Foreground,  derived  from  H  —  K  color 

This  work  . 

-21.4 

-22.0 

Foreground,  derived  from  H  recombination  lines 
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TABLE  5 

Dynamical  Mass  Estimates 


Radius  Mass 


Authors  (pc)  (109  M0)  Technique 


Shen  &  Lo  1995 .  390  2.3  CO  (1  — 0)  measurements 

Achtermann  &  Lacy  1995 .  344  1.0  [Ne  n]  12.8  /im  measurements 

Gotz  et  al.  1990  .  390  1.5  Na  D,  Ha,  H  i  21  cm 

344  1.2  Na  D,  Ha,  H  I  21  cm 


estimate  the  mass  enclosed  in  a  30"  aperture  (r  =  225  pc),  we 
use  the  mass  model  proposed  by  Gotz  et  al.  (1990), 

M(<r)  =  5  x  10V-33  +  4  x  >  (2) 

where  M(<r)  is  the  mass  (M0)  interior  to  the  radius  r  (pc). 
Using  this  model,  the  mass  within  a  30"  aperture  is  7  x  108 
M0.  This  value  is  the  same  as  that  assumed  by  Rieke  et  al. 
(1993). 

The  molecular  gas  mass  has  been  estimated  using  a 
variety  of  methods.  These  methods  are  summarized  by  Wild 
et  al.  (1992).  A  selection  of  these  measurements,  demonstrat¬ 
ing  the  large  discrepancies  in  the  derived  molecular  masses, 
is  shown  in  Table  6.  It  has  been  shown  that  the  molecular 
gas  conversion  factor  appropriate  for  M82  is  much  smaller 
than  the  galactic  value  (e.g.,  Lo  et  al.  1987;  Paper  I).  In 
addition,  we  have  also  shown  in  Paper  I  that  the  geometry 
of  the  starburst  region  can  be  approximated  by  a  molecular 
torus  surrounding  a  central  cavity  containing  predomi¬ 
nantly  ionized  gas  and  stars.  Therefore,  the  estimates  for  the 
molecular  mass  in  Table  6  should  be  much  larger  than  the 
true  value  appropriate  for  the  starburst  region.  Thus,  we 
assume  that  the  central  30"  does  not  contain  more  than  108 
Me  of  molecular  gas,  and  most  probably  contains  a  negligi¬ 
ble  amount. 

The  above-mentioned  mass  constraints  and  extinction- 
corrected  K-band  luminosity  can  now  be  applied  to  our 
stellar  population  synthesis  models.  Using  our  extinction- 
corrected  absolute  K  magnitude  of  —  22.0,  we  find  that  only 

2.5  x  10s  M0,  36%  of  the  total  dynamical  mass,  is  needed 
in  the  burst  in  order  to  accommodate  a  Salpeter  IMF  that 
extends  from  0.1  to  100  MQ,  and  that  the  remaining 

3.5  x  108  M0,  or  ~  50%  of  the  total  dynamical  mass,  can 
be  locked  up  in  an  old  population,  assuming  that  the  burst 
age  is  ~  107  yr.  The  simple  modification  to  the  two-screen 
dust  geometry  (see  §  4.4)  gives  approximately  the  same 
results.  In  order  to  illustrate  the  relationship  between  MK, 
MTot,  and  mL,  we  have  taken  the  range  of  extinction- 
corrected  K-band  magnitudes  from  Table  4  and  calculated 
the  lower  mass  limit  required  to  produce  these  magnitudes 
for  a  burst  age  of  ~  107  yr.  Figure  11  is  a  graph  of  mL  versus 
the  total  mass  assumed  in  the  burst,  MXot.  As  can  be  seen,  a 


moderate  variation  in  the  mass  allocated  to  either  the  burst, 
shown  above  to  be  an  uncertain  quantity,  or  the  extinction- 
corrected  K  flux  density — variations  that  can  be  tolerated 
by  the  current  observational  constraints — can  result  in  a 
lower  mass  cutoff  for  the  starburst  region  that  is  much  more 
like  the  solar  neighborhood  than  has  been  assumed  pre¬ 
viously.  Note  that  it  has  been  argued  that  the  actual  solar 
neighborhood  IMF  is  steeper  than  the  Salpeter  IMF  at  high 
masses  (Massey  et  al.  1995b).  However,  our  assumed  IMF 
does  not  differ  greatly  from  other  assumed  IMFs  at  high 
masses.  If  the  more  recent  Miller  +  Scalo  IMF  (Miller  & 
Scalo  1979)  is  assumed  instead  of  the  Salpeter  IMF,  nor¬ 
malizing  to  the  same  total  mass  and  assuming  the  same 
lower  and  upper  mass  cutoffs,  the  fraction  of  the  total  mass 
in  stars  with  masses  greater  than  10  Mq  is  10%  compared 
with  12%  for  the  Salpeter  IMF.  Thus,  a  Miller  &  Scalo 


Fig.  11. — Relationship  between  the  lower  mass  cutoff  and  the  total 
mass  assumed  in  the  burst  for  various  values  of  the  extinction-corrected  K 
magnitudes.  The  total  dynamical  mass  for  the  starburst  region  (450  pc)  is 
7x10  8M0. 


TABLE  6 


Molecular  Mass  Estimates 


Authors 

Radius 

(pc) 

Mass 
(108  Mq) 

Technique 

Wild  et  al.  1992  . 

800 

1.8 

ClsO  (2-1) 

Carlstrom  1989 . 

375 

1.0 

CO  (1—0),  assumed  optically  thin 

Lo  et  al.  1987  . 

700  x  200 

0.6 

CO  (1—0),  assumed  optically  thin 

Olofsson  &  Rydbeck  1984  . 

500 

0.6 

CO  (1—0),  assumed  optically  thin 

Young  &  Scoville  1984 . 

500 

4.0 

CO  (1—0),  assumed  optically  thick 

Sutton,  Phillips,  &  Masson  1983 . 

250 

0.3 

CO  (2  —  0),  assumed  optically  thin 
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IMF  extending  from  0.1  to  100  M0  can  replicate  the 
observed  extinction-corrected  K-band  magnitude  if  the 
total  mass  assumed  in  the  burst  is  increased  slightly. 

If  a  large  percentage  of  the  total  mass  is  in  the  form  of  an 
old  underlying  stellar  population,  the  lower  mass  cutoff 
required  of  the  burst  population  becomes  greater  than  1 
M0.  For  example,  if  we  assume  that  ~  70%  of  the  mass  is  in 
an  old  population,  the  lower  mass  cutoff  required  to  repli¬ 
cate  our  observed  extinction-corrected  K-band  magnitude 
is  1  Mq.  Since  it  is  unknown  how  much  mass  is  in  the  form 
of  the  underlying  stellar  population,  the  possibility  that 
almost  all  of  the  mass  is  in  this  form,  and  therefore  that  the 
burst  population  is  deficient  in  low-mass  stars,  cannot  be 
ruled  out.  Our  model  predicts  a  supernova  rate  of  0.08-0.1 
yr-1  for  ages  between  3  x  106  and  1  x  107  yr  when  our 
extinction-corrected  K  magnitude  of  —22.0  is  matched. 
This  is  consistent  with  estimates  derived  from  radio  obser¬ 
vations  (e.g.,  0.1  yr-1,  Bartel  et  al.  1987;  0.05  yr-1,  Muxlow 
et  al.  1994).  Again,  the  reader  is  reminded  that  the  purpose 
of  this  discussion  is  to  illustrate  primarily  the  relationship 
between  the  various  observed  quantities  and  the  conclu¬ 
sions  regarding  the  lower  mass  cutoff  of  the  IMF.  In  reality, 
the  treatment  of  the  integrated  properties  from  the  entire 
500  pc  as  a  single  instantaneous  starburst  is,  as  emphasized 
in  this  paper,  an  oversimplification. 

In  their  most  recent  work  on  M82,  Rieke  et  al.  (1993), 
using  their  updated  extinction-corrected  K  luminosity 
(McLeod  et  al.  1993;  see  Table  4),  have  modeled  the  stellar 
population  in  a  30"  aperture  assuming  two  Gaussian  star- 
bursts,  the  first  to  provide  the  red  supergiants  and  super¬ 
novae  for  the  X-ray  wind  and  the  second  to  produce  the  UY 
flux  (ionizing  photon  flux).  Although  their  revised 
extinction-corrected  K  luminosity  is  substantially  fainter 
than  their  1980  result  (Rieke  et  al.  1980),  their  conclusion  of 
a  low-mass  deficient  IMF  remains  unchanged.  The  crucial 
factor  behind  this  conclusion  is  the  fact  that  in  their  model, 
they  find  that  at  no  time  do  the  K  and  UV  luminosities 
match  their  observed  values.  At  early  times,  when  the  stellar 
population  is  dominated  by  hot  young  stars  and  the  model 
UY  flux  is  similar  to  that  observed,  the  K  flux  density  falls 
short  of  its  observed  value.  At  later  times,  when  more  stars 
have  evolved  into  supergiants  and  the  K  flux  density  is 
replicated,  the  UV  flux  falls  short.  Despite  their  incorpor¬ 
ation  of  two  starbursts,  either  the  UV  flux  or  the  K  flux 
density  falls  short  by  a  factor  of  2  if  a  solar  neighborhood 
IMF  is  assumed  and  the  total  mass  assumed  in  the  burst  is 
fixed  at  2.5  x  108  M0.  In  order  to  match  synchronously 
these  observed  quantities  without  increasing  the  total  mass 
assumed  to  be  participating  in  the  burst,  the  IMF  needs  to 
suppress  the  formation  of  low-mass  stars.  Part  of  the  com¬ 
plexity  of  these  models  arises  from  the  variation  in  the 
stellar  population  within  the  450  pc  region  being  modeled. 
As  seen  from  our  Bry  and  K  images  (Paper  I),  the  UV  flux  is 
high  where  the  K  flux  density  is  low  and  vice  versa.  In  fact, 
in  a  2"  aperture  centered  on  the  nucleus,  the  ratio  of  the  K 
flux  density  to  the  Bry  flux  is  a  factor  of  3  times  larger  than 
the  ratio  in  a  30"  aperture.  However,  the  ratio  of  the  dered- 
dened  K  luminosity  to  dynamical  mass  (using  eq.  [2]) 
remains  the  same  in  these  two  apertures.  A  comparison  of 
the  more  reliable  mass  estimates  derived  from  near-infrared 
absorption  features  with  those  derived  from  the  ionized  gas 
emission  lines  shows  that  LK/M  is  even  smaller  in  the 
central  few  arcseconds  of  M82  than  in  the  starburst  disk 
(Gaffney  et  al.  1993).  Thus,  a  solar  neighborhood  IMF  that 


includes  low-mass  stars  could  be  completely  consistent  with 
observations  if  the  region  modeled  were  smaller  in  extent. 

7.  SUMMARY  AND  CONCLUSIONS 

In  this  paper,  we  have  compared  a  starburst  stellar  popu¬ 
lation  model  with  a  series  of  detailed  high  spatial  resolution 
imaging  observations  described  in  depth  in  Paper  I.  Unlike 
most  previous  studies  of  M82,  our  data  have  allowed  us  to 
identify  and  model  the  properties  of  the  individual  stellar 
clusters  in  the  central  500  pc.  In  addition  to  our  published 
near-infrared  hydrogen  recombination  line  imaging,  near- 
infrared  broadband  imaging,  and  3.29  jim  dust  feature 
imaging  observations  (Paper  I),  we  have  also  carried  out 
imaging  observations  in  the  2.36  n m  CO  band  head.  Our 
major  results  are  summarized  as  follows : 

1.  The  CO  2.36  /zm  absorption  feature  is  seen  to  be  anti¬ 
correlated  with  the  Bry  emission.  The  CO  index  is  found  to 
be  0.22  at  the  nucleus  of  M82,  higher  than  that  of  normal 
galaxies  (F78).  Smaller  values  are  seen  farther  from  the 
nucleus,  suggesting,  along  with  other  star  formation  diag¬ 
nostics,  that  the  nucleus  contains  an  older  population  of 
stars  and  that  the  starburst  is  propagating  outward. 

2.  A  comparison  of  our  model  quantities  with  obser¬ 
vations  allows  us  to  determine  the  variation  in  age  within 
the  starburst  complex.  Our  analysis  implies  that  the  typical 
age  for  the  starburst  clusters  is  107  yr.  The  inferred  age 
dispersion  between  the  stellar  clusters  in  the  central  500  pc 
is  found  to  be  6  x  106  yr.  If  the  starburst  in  M82  is  propa¬ 
gating  outward  from  the  center,  this  corresponds  to  a  veloc¬ 
ity  of  propagation,  originating  in  the  center,  of  ~  50  km  s” 1. 

3.  An  analysis  of  the  integrated  properties  in  a  30"  aper¬ 
ture  (r  =  225  pc)  reveals  that  our  extinction-corrected  K 
magnitude  of  —22.0  can  be  accommodated  by  a  Salpeter 
IMF  extending  between  0.1  and  100  M0  if  we  assume  that 
36%  of  the  total  dynamical  mass  within  the  central  450  pc  is 
in  the  burst  population,  leaving  more  than  50%  of  the  mass 
for  a  preexisting  population.  Our  high-resolution  data 
reveal  that  when  the  region  modeled  is  decreased  in  extent, 
the  higher  mass-to-K  luminosity  ratio  and  the  smaller  UV 
photon  flux  make  it  increasingly  easier  to  match  synchro¬ 
nously  the  observed  properties  without  invoking  a  high- 
mass-biased  IMF. 
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